Domain wall suppression in trapped mixtures of Bose-Einstein condensates 
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The ground state energy of a binary mixture of Bose-Einstein condensates can be estimated for 
large atomic samples by making use of suitably regularized Thomas-Fermi density profiles. By 
exploiting a variational method on the trial densities the energy can be computed by explicitly 
taking into account the normalization condition. This yields analytical results and provides the 
basis for further improvement of the approximation. As a case study, we consider a binary mixture 
of 87 Rb atoms in two different hyperfine states in a double well potential and discuss the energy 
crossing between density profiles with different numbers of domain walls, as the number of particles 
and the inter-species interaction vary. 

PACS numbers: 67.85.Hj, 67.85.Bc, 03.75.Mn 



I. INTRODUCTION 

Binary mixtures of Bose-Einstein condensates [T] are 
of great interest due to their complex dynamical features 
and their role in the emergence of macroscopic quantum 
phenomena [2j. Mixtures are experimentally available 
and usually made up of two alkali atomic species [3HS]- 
They generally display repulsive self-interaction and are 
confined by various external potentials. Depending on 
the inter-species interaction, two classes of stable config- 
urations are possible: mixed and separated. The latter 
are more interesting, since it is in this case that the ob- 
servation of phenomena such as symmetry breaking and 
macroscopic quantum tunnelling of one species through 
the other one [T0HI2] is possible. Evidence of phase sep- 
aration has been observed in [9] . Many recent articles 
are devoted to the investigation of dynamical effects in 
mixtures, such as vortices and solitons (see [TMTrj] for 
recent experimental and theoretical studies). 

Different approaches are possible in order to study the 
ground state of these systems. If the number of par- 
ticles in the condensate is very large compared to the 
number of particles in the excited states, the fields asso- 
ciated to the two species can be treated as classical wave 
functions. This approach leads to the Gross-Pitaevskij 
equations |17) . which are nonlinear Schrodinger equa- 
tions obtained by finding stationary points of the zero- 
temperature grand-canonical energy of the system. The 
Thomas-Fermi (TF) approximation, which consists in 
this case in neglecting the kinetic energy of the system, is 
then usually applied [17] ■ A great deal of results are ob- 
tained in particular cases, such as confinement by a hard 
wall trap [T5] , harmonic or lattice potentials [TH] and ax- 
isymmetric traps |20j , also in the presence of the gravita- 
tional force [21]. The problem of the stability of mixtures 
has been tackled also with renormalization group tech- 
niques [2"2] . 

The presence of the kinetic term in the energy func- 
tional of a binary mixtures leads in particular to the reg- 
ularization of possible domain walls, which sharply sepa- 



rate the two species in the TF ground states [T01 . 
This is generally related to the problem of minimizing 
the surface energy, that is also found in the theory of 
superconductivity [27] . 

In this work we introduce a variational method in order 
to approximate the Gross-Pitaevskij solution in a neigh- 
borhood of a domain wall and estimate the total energy 
of a mixture. Our technique explicitly takes into account 
the normalization of the condensate wave functions and 
ensures complete analytical feasibility. 

This article is organized as follows: in Section [Tl] we 
summarize results obtained in the TF approximation 
which are relevant to our analysis; in Section |III| we in- 
troduce the regularization technique and obtain results 
regarding the energy increase with respect to the TF ap- 



proximation; in Section IV we examine a case study in 
which macroscopic effects related to domain wall suppres- 
sion can be observed and present a quantitative phase 
diagram; in Section [V] we suggest further possible devel- 
opments of our technique. 



II. BINARY MIXTURES 



A. Energy functional 



Let us consider a binary mixture of Bose-Einstein con- 
densates, confined by the external potentials Vk(x), with 
k = 1,2. We assume that the particles are tightly con- 
fined in the transverse directions, so that the system is 
quasi one-dimensional |17j . Let U^k > be the parame- 
ters that determine the self interaction between particles 
of each species and Uu = U21 > the inter-species in- 
teraction parameter. The ground state of the system 
is determined by the coupled Gross-Pitaevskij equations 
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with k = 1,2. They are the variational equations of the 
quartic energy functional 
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Here ipk(z) are the condensate wave functions, whose 
squared moduli p k (x) = ^(x)! 2 represent the local den- 
sities of each species, N k are the numbers of particles 
making up the condensates, while 71- , V k and U are the ki- 
netic, potential and interaction energy, respectively. We 
will assume henceforth that the condensate wave func- 
tions are real, since energy minimization requires their 
phases to be constant. 



B. Thomas- Fermi approximation 

The TF approximation |17j is in this case equivalent 
to neglecting the kinetic terms in Q-Q, so that 

£tf(pi,P2) = ^2 v fe (7^:) +u( A /p^,^). (5) 



The value of the adimcnsional parameter 



(6) 
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is crucial in qualitatively determining the TF ground 
state. If a > 1, which is the case of interest here, the 
ground state density profiles are completely separated in 
adjacent regions divided by domain walls. In the regions 
where only species k is present, the solution to the TF 
equations reads 
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U, 



A- A- 



(7) 



Equation ^ completely determines the functional de- 
pendence of the densities on the external potential (and 
on the chemical potentials), once their supports are 
given. Here we are interested in the ground state solu- 
tion, so that the supports are determined by minimizing 
the energy of the system ^ as a function of the number 
and the positions of the domain walls. In Ref. [55] we 
proved that for continuous trapping potentials stationar- 
ity of energy requires that the densities at a domain wall 
at Rj satisfy 



2Pl{Rj)- 



(8) 



Moreover, in the special case V(x) = V\{x) — V 2 {x), 
conditions ([8| imply that the external potential, and thus 
the density of each species, should be the same at all 
domain walls. 

The TF approximation works very well for large num- 
bers of particles, and provides a good estimate of the en- 
ergy of the system. Despite being small with respect to 
the potential energy, corrections due to the kinetic term 
give nonetheless rise to macroscopic effects. The most 
relevant of such effects is the crossing between station- 
ary states with different numbers of domain walls in the 
ground state [TTJ [12] . It is thus necessary, in order to de- 
termine the actual ground state of a mixture, to consider 
a regularization scheme of the TF density profiles, that 
enables to smooth parts of the TF profiles, like domain 
walls and zeros, that provide an infinite contribution to 
the kinetic energy [17) . 



III. VARIATIONAL REGULARIZATION OF 
DENSITY PROFILES 

Many attempts have been made so far in order to con- 
sistently estimate the energy corrections due to the ki- 
netic contribution. The pioneering works by Ao and Chui 
[10] and Timmermans [23] are based on an exponential 
approximation of the TF density profiles, by extrapo- 
lating the solution of the Gross-Pitaevskij equations far 
from a domain wall, and on a linear approximation of the 
regularized walls, respectively. Other authors 24, 26J put 
forward a regularization with fixed chemical potential in 
various regimes. 

The task we will try to accomplish in this work is to 
find a proper domain wall regularization, which provides 
a reliable approximation to the ground state profile and 
energy of a binary mixture in a trapping potential, and 
which is at the same time an analytically manageable 
trial function. The approximations in [ 101 [23] are based 
on rather crude ansatz, but provide a good estimation of 
the order of magnitude of the energy changes due to the 
kinetic terms. However, we will try to find better upper 
bounds to the ground-state energies, since it is possible 
that a small energy change results in macroscopic dif- 
ferences in density profiles. Our approximation will be 
referred to a system with fixed numbers of particles, and 
will strictly rely on the preservation of the normalization 
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FIG. 1: (Color online) Plot of the trial density profiles |9f- 
(fOl in a neighborhood of a domain wall placed at Ro- The 
solid (red) line represents the density of the first species, the 
dashed (black) line the density of the second species, with 
Uu < t/22. 



conditions. It is indeed difficult to use results obtained 
with fixed chemical potential [2H US] in this case, since 
for non trivial external potentials it is generally impossi- 
ble to invert the normalization conditions (|3| and explic- 
itly express the chemical potentials as functions of the 
numbers of particles. 



A. Domain walls 

In order to regularize the domain walls, it is sufficient 
to replace the singular TF density profiles with a contin- 
uous function, with bounded first derivative. The mini- 
mization of the kinetic energy leads to tails of each species 
penetrating through the domain wall. 

Our choice of trial profiles is based on an exponential 
tail regularization of the TF solutions. Let us consider a 
domain wall placed at position i? - In its neighborhood 
we assume as profiles the continuous functions 



-2(a;-.Ri)/Ai 



if 
if 



x < Ri 
x>Ri 



and 



P2{x) 



,TFm_\ 2(x-R 2 )/A 2 



if 
if 



x > R 2 
x<R 2 



(9) 



(10) 



with Ri < R < R 2 , A fe > 0, and pJ F the TF densi- 
ties ([7]). See Fig. [I] The points Rk are fixed in order 
to locally preserve the normalization conditions. This 
is accomplished by imposing that the integral of the re- 
moved part of the TF profiles be equal to that of the new 
exponential tails: 
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which gives 
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In the following, we will neglect all corrections depend- 
ing on the first derivatives of the TF densities, p TF ' = 
dp TF /dx, which depend linearly on the derivatives of the 
external potentials. This assumption, which involves con- 
ditions on the (by now arbitrary) parameters A/j, will 
be expressed at the end of our calculations in term of 
physical quantities. It is remarkable that under this ap- 
proximation the regularized density profiles still satisfy 
condition ^ at Rq. 

Once the trial profiles are chosen, we proceed to the 
computation of the difference in potential and interaction 
energy with respect to the TF densities. It is convenient 
to explicitly show the dependence of the energy on the 
largest penetration length, say A 1; and the ratio r\ = 
^-2/^-1 < 1- Since a > 1, the mixing must result in an 
increase in potential energy. The self interaction energy 
associated to the exponential tail of p\{x) reads 



h(x) dx = -Atipi* (Ro)) 



(13) 



This contribution replaces the self interaction energy of 
the removed TF density: 
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The same results hold for the second species after the 
substitution Ai — > 77 Ai. As expected, extending the den- 
sity profiles implies a reduction in the self interaction 
energy, U sc \f =U\\ + U 22l which reads 
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A positive contribution comes from the inter-species in- 
teraction terms, Ui ntcr = U\ 2 + which are due to the 
penetration of the tails in the bulk of the other species 
and their superposition around i?o- (Remember that in a 
TF separated configuration Wi ntC r = 0-) The total change 
in the inter-species interaction energy reads 
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Observe that the limit 77 — > 1 is finite. It is easy to verify 
that no corrections come from the interaction with the 
external potential at the chosen order of approximation. 



The results ( 15 )-( 16 ) are found under the hypothesis that 



the distance separating the considered domain wall from 
other possible walls is much larger than the A^'s. 

Let us now consider the contributions from the kinetic 
energy. The value of the kinetic energy of the densities 
in the bulk is consistently neglected in our approxima- 
tion, being 0((k k p TF ' / p TF ) 2 ) with respect to the lead- 
ing terms. On the other hand, the contribution across the 
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domain wall depends on the inverse penetration lengths 
and reads 
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1 + ^ pf(iio), (17) 



is the ratio between the healing lengths 
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(A; = 1,2) of uniform condensates whose densities are 

The energetic contributions ( 15 )— ( 17 1 depend on the 
free parameters A k . Heuristically, as a benchmark, one 
can evaluate them at A k = \J a — 1, the penetration 
lengths associated with the exponential tail of one species 
into the bulk of the other one, which are obtained via 
an approximated Gross-Pitaevskij equation [TO]. Besides 
yielding a divergence of the penetration for a close to 1, 
which is unphysical for a system with finite numbers of 
particles, this is not the optimal choice, and it can even 
yield an energy 20% larger than the minimum for physical 
parameters. In a situation where small energy changes 
are involved, this discrepancy can be very relevant. The 
best strategy is to minimize the total energy of the wall 
over the family of trial functions (|9|-(10) parametrized 
by Ai and 77 = A2/A1: 
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This will enable us to get a much more accurate upper 
bound for the ground state energy of the binary mixture. 
The minimum in ( 20 1 is attained at a single point 



(Ai (a), 77(a)) and reads 
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The optimal ratio between the penetration lengths, 77(a), 
is the solution to the transcendental equation 
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FIG. 2: (Color online) Behavior of the function $ a (77(a)), 
obtained by solving Eq. ( |23[ ), for different values of 770- Solid 
(blue) line: 770 = 0.99; dashed (red) line: 770 = 0.75; dotted 
(green) line: 770 = 0.5. 
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FIG. 3: (Color online) Comparison between the min imum 
energy [units of £,iUn(pi F (Rq)) 2 ] of a domain wall |2l| (solid 
blue line) and the energy computed for penetration lengths 
£,k/V a — 1> f° r 770 = 0.9 (dashed red line). 



In Fig.[2]the function $ Q (77(a)) is plotted versus a for dif- 
ferent values of the parameter 770, while in Fig. [3] a com- 
parison is displayed between the minimized energy and 
the result computed with the bulk penetration lengths 
£fc/Vo--T. It can be observed that the minimizing en- 
ergy has a monotonic behavior with a. The optimal 
penetration length 
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17(a) 



1 (r?(a)) ' 



(24) 



is plotted versus a (for 770 = 0.9) in Fig.|4j The conditions 
ensuring that terms depending on the first derivatives of 
the densities can be neglected are thus summarized by 
the following inequality 
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the form tanh(a;/-\/2£fc) and to a saturation of the do- 
main wall energy. Since the typical variation lengths of 
the densities are in this case V^Cfc, we can bind the valid- 
ity of our ansatz to values of a verifying A&(a) < V%£k 
(see Fig. [4] for a comparison) , typically corresponding to 
values a < 2.5 -j- 3, which matches well the experimental 
ranges. 

B. Profile edges 
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FIG. 4: (Color online) The solid (blue) line represents the 
optimal penetration length Ai(a) (units of £i) plotted against 
a for r/o = 0.9. The dashed (red) line is the bulk penetration 
length £i/\/a — 1 |10| . plotted here for comparison. 



domain walls is always larger than Afc's, thus validating 
the results in Eqs. (p~5]) - (fT6l) . 



Since the dependence on a and p{ in (21) is factor- 



ized, the total correction for a TF configuration with do- 
main walls placed at {Ri, . ■ . , R w } is 

w w 

AfW =J2^n(Rj) = C^ipTiR,)?' 2 , (26) 



with 



C a = ^(i7 11 /2m 1 ) 1 / 2 $ Q (??(a)). 



(27) 



For confining potentials proportional to each other, in- 
cluding the interesting case in which they are equal, the 
energy increase (21) is independent of the specific do- 
main wall, since ( 8 ) implies that the densities are the 
same at all the walls of a stationary configuration, namely 
pJ F (Rj) = pJ F (R ) for any j [29]. This implies that the 
energy correction to a TF configuration with w domain 
walls is simply 



ASM = w£ wall (R ) 



On the other hand, the condition of applicability of 



our approximations (25) depends on the specific domain 



wall through the first derivative of the potentials at each 
wall. Of course, a uniform control on the derivatives of 
the potentials would give a sufficient condition for their 
applicability to all possible configurations. 

It shoud be emphasized that the validity of our ap- 
proximation has an upper bound in a. Indeed, it can 
be deduced from Eq. ( [22] ) and Fig. [3] that the function 
$ Q (77(a)) is not bounded from above as a increases, lead- 
ing to a divergent correction to the TF energy, which 
frustrates regularization attempts. The reason lies in the 
fact that for a — > 00 the exponential tail ansatz is no 
longer justified, since in this case one species feels the 
other one like a hard wall, thus leading to solutions of 



The presence of first order zeros at the edges of the TF 
density profiles, which are compactly supported, gives 
rise to logarithmic divergencies in the kinetic energy. 
Even this situation can be tackled by a proper regu- 
larization of the densities, leading to an increase in the 
potential energy and to a kinetic contribution |17j . To 
this end, we consider a zero at position xq and con- 
ventionally consider a TF profile pJ F (x)9(xq — x) in a 
neighborhood of Xq, 6 being the unit step function. The 
regularization is based on the solution of the linearized 
single-condensate wave function in a neighborhood of xq, 
namely y/~p~k oc Ai(x/S), with Ai(y) the proper Airy func- 
tion [30J , decreasing as y — > 00, and 



6k = 



( 2m k F k (x ) 



h 2 
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dV k 
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(29) 

the characteristic length. Corrections to a solution based 
on the linear approximation of the potential are negligible 
if the second derivative of the potential is much smaller 
than {2m k F^(x )/h 2 ) 1/3 - 

Following such scheme, we consider the trial family of 
(continuous and positive) regularized TF profiles 



Pk(x) 



if x < x k 

2 . f . , (30) 
it x > x k 



with f(y) = Ai(y)/Ai(0). As in the case of the domain 
wall, the point x k is determined by requiring the local 
normalization condition to be fulfilled. The result, anal- 



(28) ogous to (12), reads 



Xk - Xq = -2I $k, 



(31) 



where Iq = J °° dyf 2 (y) ~ 0.53. Summing the variations 
of the self-interaction energy (negative) and of the ex- 
ternal potential energy (positive) to the kinetic energy 
yields the total energy change due to the regularization 
of the zero of the density profile 



e£L(x ) ~ 0.274 



m k Ukk 



Fk(x ). 



(32) 



This contribution depends on the specific zero, as well as 
on the species k, since the first derivatives of the potential 
at its zeros are generally not related. 

In order to get a feeling for the orders of magnitude of 
the various energies, one can consider the simple case in 
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which the trapping potential is well approximated by a 
power law 



V(x) 



(33) 



A TF zero, placed at xq, of the density profile of species 
k, is determined by the condition /i^ = V(xq) ~ Xq. 
By considering the last relation and the normalization 
conditions, one can obtain the scaling law of the chem- 
ical potential with respect to the number of particles, 
namely p k 

ing laws of the TF energy, £tf 



one can obtain the seal- 

^(an+l)/(n+l) t and 



of the kinetic energy, including the contribution ( 32 ) , 
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energy of a domain wall de- 



pends on (p 



TF^3/2 
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and thus £ wa ii ~ N, 



3n/(2n+2) 



The higher order corrections of 0(Akp' k / pk) terms acci- 
dentally scale like the kinetic energy. Thus, in a config- 
uration with w domain walls, the trial ground state has 
energy 



£ = Utf + w£ wa u + O 



where 

is the energy of the TF densities. 



(34) 



(35) 



IV. DOMAIN WALL SUPPRESSION: A CASE 
STUDY 

As an application of the previous results we consider 
in this section the energy crossing between configurations 
with different numbers of domain walls in a double well 
potential. We will consider a physical situation in which 
a crossing between ground states with a maximal and a 
minimal number of domain walls can be observed. We in- 
troduce an operational way to control the crossing, based 
on a scaling property of the TF energy functional, thus 
suggesting a possible experimental realization. 

Notice first that the effect of the kinetic energy on the 
ground state of a mixture in a square well potential is 
trivial, since in this case the TF energy for separated 
configurations depends only on the volumes occupied by 
the two species and not on how they are distributed inside 
the well. Therefore, inclusion of the domain wall energy 
immediately enables us to identify the configuration with 
a single domain wall as the ground state. 

The situation is much more interesting for potentials 
that vary over the region occupied by the mixture. In 
the TF theory, density profiles with a maximal number 
of domain walls are usually energetically favored, espe- 
cially when the ratio of the self-interaction coefficients 
is very close to one [29] , However, the inclusion of the 
kinetic energy can drastically change this picture. Each 



domain wall has an energetic cost, expressed by (21 1, 



whose effect on the total energy decreases as the num- 
bers of particles increase. Thus, for very large number 
of particles, the ground state is more likely to have a 
maximal number of domain walls, but if the number of 
particles decreases or the a parameter increases, it can 
become more convenient to reduce the number of walls. 

As a test ground for the effectiveness of our method, 
let us consider an example of this phenomenon, that was 
analyzed by numerical integration of the coupled Gross- 
Pitaevskii stationary equations in Ref. [TT]. The poten- 
tials are harmonic and the possible competing ground 
states have one or two walls: a mixture of 87 Rb atoms 
with mi = 77i2 = m = 1.45 x 10 - 25 kg in two hypcrfine 
states was considered, with scattering lengths a% = 5.36 
nm and a 2 = 5.66 nm. The longitudinal trapping fre- 
quency was fixed to uj = 2-7T x 90 Hz, with the transverse 
trapping frequency 30 times larger. The authors were 
able to build a phase diagram showing the crossing be- 
tween single- and double-wall configurations by applying 
numerical techniques. An explicit comparison of the re- 
sults for the total energies of the configurations are given 
in [UJ for Nx = N 2 = 2000 and a = 1.18. 

In Table [IJ we compare our results with those of 
Ref. [TT]. It is manifest that, while for the symmetry- 
preserving (double-wall) state the two results are iden- 
tical up to the fourth significant digit, our regulariza- 
tion method enables us to attain a stricter upper bound 
for the ground-state energy of the symmetry-breaking 
(single-wall) configuration. Since the choice of the trial 
densities in the energy functional is based on the physics 
of the phenomenon, it is not surprising that our ana- 
lytical regularization technique, together with the ex- 
act results coming from TF, leads to a better approx- 
imation of the ground state of a binary mixture than 
the accurate numerical integration of the coupled Gross- 
Pitaevskii equations [TT] , 

The approximation on the energy of the trial densi- 
ties was proved to be robust by a numerical check, in 
which the total energy of the regularized TF profiles is 
computed by numerical integration, showing only a slight 
increase in the fifth significant digit. The differences in 
the estimate of the ground-state energy of the symme- 
try breaking configurations, which are general and not 
restricted to the aforementioned case, lead to a different 
phase diagram in the plane (Ni = N 2 = N,a), shown in 
Fig. [5] that should be compared with that in Fig. 2 of 
Ref. [TT] . In our case the transition line is shifted by a 
factor ~ 1.5 with respect to the N axis (towards larger 
values of N). Thus, according to our analysis, the sym- 
metry breaking ground state is present in a larger region 
of the (N, a) plane, where it was previously not expected. 



A. Scaling properties 

It is convenient to study the scaling properties of the 
energy terms in ^ under a dilation. If lengths scale 
as x — > xja with a > 0, it is easy to see from Q-Q 
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TABLE I: Total energy of a binary mixture of Rb atoms in a harmonic potential with longitudinal frequency lj = 2-k x 90 
Hz and transverse frequency 30 times larger, for N 1 = N% = 2000 and a — 1.18. The results obtained with the analytical 
approximation schemes proposed in this article are compared with those numerically obtained in Ref. [11] , 

One wall (symmetry breaking) Two walls (symmetry preserving) 
Wtf 87.091 86/772 

£ = Wtf + A£ (ro) 87.486 87.423 

£ numerically computed in [11] 87.551 87.426 
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FIG. 5: (Color online) Ground state phase diagram in the N— 
a plane for a binary mixture of 87 Rb atoms in a harmonic po- 
tential with longitudinal frequency w = 90Hz and transverse 
frequency 30 times larger. In the light grey region (yellow in 
the online version) the ground state is a symmetry preserving 
configuration with two domain walls, while in the dark grey 
region (red in the online version) a symmetry breaking con- 
figuration, with a single domain wall, is energetically favored. 



unions of intervals whose ends are domain walls or edges, 
located at {Rj} = V' 1 ^}) and {Pj fc) } = V -1 (-{>*,})> 
respectively. Observe now that if the potential is scaled 
as 

W(x) = V(x/a), (37) 

while leaving unchanged the chemical potentials jji k 
and the domain wall potential v, the density pro- 
files a k VF (x) = pJ F (x/a) are still TF solutions, cor- 
responding to the potential W, to supports afl k = 
Qk({aRj}, {aP^}), and to numbers of particles 

f dyal v {y) =af dx P J F (x) = aN k . (38) 

Moreover, the energy of the scaled configuration is re- 
lated to the previous one by 

U TF {aN k ,W,w)=aU T F (N k ,V,w). (39) 

On the other hand, the energy contribution of the do- 
main walls is unchanged by the scaling, since it depends 
only on fixed quantities, namely the number of walls, the 
chemical potentials and the domain wall potential. 

B. Energy crossing 



that the kinetic energies scale as T k — > T k /a, while the 
potential and interaction energies scale as V k — > aV k and 
hi — > aU, and accordingly the numbers of particles scale 
as Nk — >• aN k . Therefore, the larger a and the numbers 
of particles, the smaller the ratio between kinetic and 
potential energy. Let us look at this property in more 
details. 

Consider TF density profiles p k F in a separated con- 
figuration given by ([7| and fix their supports by choosing 
the number w of domain walls and their positions Rj , all 
satisfying V(Rj) = v with v a constant value. Let the 
density profiles be normalized to N k . The TF potential 
energy of the separated configuration reads 

U TF (N k ,V,w) 

= W dx (v{x)pF(x) + ^pF(*f) ■ (36) 

k=l,2 Uk ^ ' 

The integration domains Q k = £l k ({Rj}, {P^}) are 



Let us consider for definiteness a physical system with 
equal numbers No of particles of the two species in 
a potential V, with w domain walls, potential energy 
Utf(No,V,w) and total domain wall correction A£^ w \ 
proportional to w. If the numbers of particles are in- 
creased to N — aNo with a > 1 and the potential is 
stretched to V(xNq/N), the TF energy of the new con- 
figuration reads 

U TF (n,v(x^ ,w\ = ^U TF (N ,V(x),w), (40) 

while the domain wall contributions remain the same. We 
can conveniently consider Nq such that the bulk kinetic 
energy and the energy corrections due to the zeros of the 
density profiles are negligible with respect to both the 
TF energy and the domain wall energy. A fortiori, they 
will be negligible for all N > No, since the bulk kinetic 
energy scales like Nq/N. 

We consider now an alternative configuration, in which 
the number of domain walls is w' . The potential energy 



Utf(N, V(xN /N),w') obeys the same scaling law (40). 
If a crossing between the total energies of the configura- 
tions exists, it occurs for a number of particles 



N*, (a)=jV - 



'U TF (N , V(x),w>) - U TF (N , V(x),w) ' 

(41) 

which is meaningful only if the differences of the potential 
energies and of the domain wall corrections have opposite 
signs. Furthermore, physical meaning can be attributed 
to the crossing only if N^ w , > Nq, since the validity of 
TF approximation is not assured for N < Nq. 

Binary mixtures of 87 Rb atoms are experimentally 
available 0(5], with mass mi = m 2 —m— 1.45 x 10~ 25 
kg, in states \F — l,mj? = +1) and \F = 2,mp = —1), 
whose s-wave scattering lengths are respectively a\ = 
5.36 nm and a 2 = 5.66 nm. The inter-species scattering 
length ai2 is tunable by approaching a Feshbach reso- 
nance [3T] (see [H 03 [S] for experimental realizations). 
Let us suppose that such a mixture is confined in a 
deformed harmonic trap, with a longitudinal frequency 
ui = 27r x 0.7 Hz, corresponding to a trapping length 
a£ — yjh/ (muje) — 1.29 x 10~ 5 m and a transverse fre- 
quency ujj_ = 500 uj£, such that aj_ — 0^/(10-^/5). Since 
we want the transverse degrees of freedom to be freezed, 
the number of particles per species N has to satisfy 
A*o aj/(a± min(ai, a 2 )) ~5x 10 4 . Moreover, in order 
to ensure the applicability of one-dimensional TF approx- 
imation, the condition Nq 3> a\/(ai max(ai, a%)) ~ 10 
must hold [17J. A good choice is then JVo = 5x 10 3 . It 
is readily verified that if the potential is (longitudinally) 
scaled as in (37), the assumption of one-dimensionality 
and the TF approximation continue to be valid. The 
one-dimensional self-interaction parameters read 

2H 2 a k 



Ukk — 



m k a\_ 



(42) 



and their ratio U11/U22 = 0,1/0,% is very close to one. 

In order to obtain a double well in the region where 
the condensates are trapped, we add to the longitudinal 
harmonic potential a cosine potential, so that 



V(x) 



~2~ 



A cos(Bx), 



(43) 



with A/(hu e ) = 107.75 and Ba e = 6.44 x 10~ 16 . 

The TF stationarity condition (jSJ), together with the 
normalization conditions, is satisfied by the seven differ- 
ent configurations represented in Fig. [6j together with 
their TF and domain wall energies. The number of do- 
main walls ranges from one to four. The configuration 
with four domain walls and with the less-self-interacting 
species placed in the minima of the external potential 
is, as expected, the minimizer of the TF energy. How- 
ever, its domain wall energy is much larger than that of 
the configuration with a single domain wall, where each 
condensate occupies on e w ell. For instance, a = 2 yields, 
according to definition ((111, ATf 4 (2) = 1.4 x 10 4 . We have 
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FIG. 6: (Color online) Lowest energy configurations in the 
double well potential (43|, for N = 5 x 10 3 and a = 2. In 



all figures the (linear) densities of species 1 (solid red line) 
and species 2 (black dashed line) are plotted vs the linear 
coordinate x. Below each figure, u = U(No,V,w)/No is the 
TF potential energy per particl e, w hile S is the specific energy 
of the domain walls as in Eq. (261, both in units hue- 



observed that for all a the only competing ground states 
are the aforementioned configurations, with one and four 
domain walls. The four-wall profile is the ground state 
only for N > N* 4(a), while for smaller the one with a 
single domain wall is energetically favored. In Fig. [7] the 
ground state phase diagram is displayed. The transition 
line is the graph of the function 



AT* 4 (a) cx & a (fj(a)) 



(44) 



and gives direct information about the function 
<!> a (fj(a)). Thus, using Eqs. plj) and p4|, the ground 



state phase diagram can be used to obtain information 
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FIG. 7: (Color online) Ground state phase diagram for a mix- 
ture with equal number of atoms Aq — N2 — N in the poten- 
tial V(x) of Eq. (431. In the light grey region (yellow in the 
online version) configurations with the maximal number of 
walls are energetically favored, while in the dark grey region 
(red in the online version) the ground state has a single do- 
main wall. The (blue) transition line represents the function 
JV?, 4 (a). 



about the domain wall energies and the optimal penetra- 
tion lengths. 



V. CONCLUSIONS AND OUTLOOK 

We discussed a variational method that yields a very 
good approximation for the total energy of a binary mix- 
ture of Bose-Einstein condensates in a separated con- 
figuration. The method is reliable and accurate. We 
have seen that, in some cases, density profiles with a 
large number of domain walls can be energetically fa- 
vored with respect to those with fewer domain walls, in 



particular when the interaction ratio a is close to one. 
At present, there is a variety of methods to find approxi- 
mate solutions to the Gross-Pitaevskii equations, ranging 
from analytical techniques |32| to numerical ones, includ- 
ing imaginary-time schemes [331 OS] and finite-difference 
methods [35U37j . It is worth emphasizing that our ap- 
proach is not aimed at solving the Gross-Pitaevskii equa- 
tions in the most general cases, but is rather optimized 
at uncovering the ground state properties of the mix- 
ture, with minimal numerical help. This enables us to 
give an immediate physical interpretation of the results, 
and makes it possible to explain in very general terms the 
observed phenomena, to predict new ones, and possibly 
to develop more refined techniques in order to extend the 
validity of our approximations. 



All results are analytical and therefore provide solid 
ground to improve the approximations. Indeed, start- 
ing from Eq. (21) and from previously obtained TF re- 
sults 



it is possible to compute, for example, correc- 
tions due to the first (finite) derivatives of the potentials, 
as well as possible domain-wall displacements and small 
variations of the chemical potentials. 

The results obtained in this paper can have practi- 
cal applications. For example, if one compares an ex- 
perimental phase diagram with the theoretical predic- 
tion, physical properties of the mixture can be estimated 
from the transition line, that depends on the ratio of the 
masses and the interaction parameters. These results can 
also help in analyzing dynamical phenomena, such as vor- 
tices and solitons, which mostly appear as perturbations 
of a stationary background, and lead to very subtle en- 
ergy changes (see e.g [TBI [38l [39]): in order to correctly 
analyze these changes, an accurate estimate of the back- 
ground energy is needed. 
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